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Abstract 

Starting from the recently obtained 2PN accurate forms of the energy and 
angular momentum fluxes from inspiralling compact binaries, we deduce the 
gravitational radiation reaction to 2PN order beyond the quadrupole approx- 
imation - 4.5PN terms in the equation of motion - using the refined balance 
method proposed by Iyer and Will. We explore critically the features of their 
construction and illustrate them by contrast to other possible variants. The 
equations of motion are valid for general binary orbits and for a class of coor- 
dinate gauges. The limiting cases of circular orbits and radial infall are also 
discussed. 

PACS numbers: 04.25.Nx, 04.30.Db, 97.60. Jd, 97.60.Lf 
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I. INTRODUCTION 



Inspiralling compact binaries are the most promising sources of gravitational radiation 
in the near future for ground-based laser interferometric detectors like LIGO VIRGO 
0, GEO600 [[| and TAMA 0. The method of matched filtering will be employed to search 
for the inspiral waveforms and extract the information they carry For this method to 

be successful, one needs to use templates that are extremely accurate in their description of 
the evolution of the orbital phase, which in turn requires a detailed understanding of how 
radiation damping (reaction) influences orbital evolution |7HlO . 



The idea of a damping force associated with an interaction that propagates with a finite 
velocity was first discussed in the context of electromagnetism by Lorentz ||11|| . He obtained 
it by a direct calculation of the total force acting on a small extended particle due to its self- 
field. The answer was incorrect by a numerical factor and the correct result was first obtained 



by Planck [12| using a 'heuristic' argument based on energy balance which prompted Lorentz 



|T3[ to re-examine his self-field calculations and confirm Planck's result, 
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where V{ is the velocity of the particle. The relativistic generalization of the radiation reaction 
by Abraham based on arguments of energy and linear momentum balance preceded by a 



few years the direct relativistic self-field calculation by Schott [T5[ and illustrates the utility 
of this heuristic, albeit less rigorous, approach |16| . 



The argument based on energy balance proceeds thus: A non-accelerated particle does 
not radiate and satisfies Newton's (conservative) equation of motion. If it is accelerated, it 
radiates, loses energy and this implies damping terms in the equation of motion. Equating 
the work done by the reactive force on the particle in a unit time interval to the negative of 
the energy radiated by the accelerated particle in that interval (Larmor's formula) the reac- 
tive acceleration is determined and one is led to the Abraham-Lorentz equation of motion 
for the charged particle. The direct method of obtaining radiation damping, on the other 
hand, is based on the evaluation of the self-force. Starting with the momentum conserva- 
tion law for the electromagnetic fields one rewrites this as Newton's equation of motion by 
decomposing the electromagnetic fields into an 'external field' and a 'self-field'. Expanding 
the self-field in terms of potentials, solving for them in terms of retarded fields and finally 
making a retardation expansion, one obtains the required equation of motion when one goes 



to the point particle limit [17 



As in the electromagnetic case, the approach to gravitational radiation damping has been 
based on the balance methods, the reaction potential or a full iteration of Einstein's equa- 



tion. The first computation in general relativity was by Einstein [|18j who derived the loss 
in energy of a spinning rod by a far-zone energy flux computation. The same was derived 
by Eddington [|19|] by a direct near-zone radiation damping approach. He also pointed out 
that the physical mechanism causing damping was the effect discussed by Laplace PD[ , that 
if gravity was not propagated instantaneously, reactive forces could result. An useful devel- 
opment was the introduction of the radiation reaction potential by Burke |^T[ and Thorne 
|22|| using the method of matched asymptotic expansions. In this approach, one derives 



the equation of motion by constructing an outgoing wave solution of Einstein's equation in 
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some convenient gauge and then matching it to the near-zone solution. Restricting atten- 
tion only to lowest order Newtonian terms and terms sensitive to the outgoing (ingoing) 
boundary conditions and neglecting all other terms, one obtains the required result. The 
first complete direct calculation a la Lorentz of the gravitational radiation reaction force was 



by Chandrasekhar and Esposito ||23|| . Chandrasekhar and collaborators p^ , p5| developed a 
systematic post-Newtonian expansion for extended perfect fluid systems and put together 
correctly the necessary elements like the Landau-Lifshiftz pseudotensor, the retarded po- 
tentials and the near-zone expansion. These works established the balance equations to 
Newtonian order, albeit for weakly self-gravitating fluid systems. The revival of interest in 
these issues following the discovery of the binary pulsar and the applicability of these very 
equations to binary systems of compact objects follows from the works of Damour |T6|j26[| 
and Damour and Deruelle Wt 



In the context of the binary pulsar timing, the accuracy reached by the Newtonian 
balance equations is amply adequate. The case of inspiralling binaries as sources for the 
interferometric gravitational wave detectors is very different. The extremely high phasing 
accuracy requirement makes mandatory the control of reactive terms way beyond the New- 
tonian. This has prompted on the one hand, work on generation aspects to compute the 
far-zone flux of energy and angular momentum carried by gravitational waves and on the 
other, work on the radiation reaction aspects to compute the effect on the orbital motion 
of the emission of gravitational radiation. As in the electromagnetic case, the computation 
of the reactive acceleration assuming balance equations is simpler than the computation 
of the damping terms by a direct near-field iteration. The computation of the energy and 
angular momentum fluxes at the lowest Newtonian order (quadrupole equation) requires the 
equation of motion at only Newtonian order. Assuming the balance equations one can infer 
the lowest order (2.5PN) radiation damping whose direct computation, as mentioned before, 
requires a 2.5PN iteration of the near-zone equations. Similarly, the computation of the 1PN 
corrections to the lowest order quadrupole luminosity requires the 1PN accurate equations 
of motion, but is potentially equivalent to the 3.5PN terms in the equation of motion. This 



motivated Iyer and Will( IW ) P8| , ^9| to propose a refinement of the text-book |30[ treat- 



ment of the energy balance method used to discuss radiation damping. This generalization 
uses both energy and angular momentum balance to deduce the radiation reaction force 
for a binary system made of nonspinning structureless particles moving on general orbits. 
Starting from the 1PN conserved dynamics of the two-body system, and the radiated energy 
and angular momentum in the gravitational waves, and taking into account the arbitrariness 
of the 'balance' upto total time derivatives, they determined the 2.5PN and 3.5PN terms in 
the equations of motion of the binary system. The part not fixed by the balance equations 
was identified with the freedom still residing in the choice of the coordinate system at that 
order. Thus, starting from the far- zone flux formulas, one deduces a formula that is suitable 
for evolving general orbits of compact binaries of arbitrary mass ratio and that includes 1PN 



corrections to the dominant Newtonian radiation reaction terms. Blanchet |3l|j32| , on the 
other hand, obtained the post-Newtonian corrections to the radiation reaction force from 
first principles using a combination of post-Minkowskian, multipolar and post-Newtonian 
schemes together with techniques of analytic continuation and asymptotic matching. By 
looking at "antisymmetric" waves - a solution of the d'Alembertian equation composed of 
retarded wave minus advanced wave, regular all over the source - and matching, one obtains 
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a radiation reaction tensor potential that generalizes the Burke-Thorne reaction potential 
33fl , in terms of explicit integrals over matter fields in the source. The validity of the balance 



equations upto 1.5PN is also proved. By specializing this potential to two-body sytems, Iyer 
and Will |29| checked that this solution indeed corresponds to a unique and consistent choice 
of coordinate system. This provides a delicate and non-trivial check on the validity of the 
1PN reaction potentials and the overall consistency of the direct methods based on iteration 
of the near-field equations and indirect methods based on energy and angular momentum 
balance. 

As emphasized earlier, much better approximations are needed to reach the precision 
of future gravitational- wave astronomy |7[]. In the limit where one mass is much smaller 
than the other, numerical and analytical computations based on black hole perturbation 
theory have been performed to the 5.5PN order |34]^39f, a recent result being the analytical 
expression to 5.5PN order for the energy flux from a test particle moving in a circular 
orbit around a Schwarzschild black hole P^fl . Ryan |^D,5T] has investigated the effect of 
gravitational radiation reaction, first on circular, and later even for non-equatorial orbits 
around a spinning black hole. Recently Mino, Sasaki and Tanaka [f|2| have derived the 
leading order correction to the equation of motion of a particle which presumably describes 
the effect of gravitational radiation reaction by two methods: an extension of the Dewitt- 
Brehme formalism and the method of asymptotic matching. 

On the other hand, for bodies of comparable masses, recently two independent teams 
43|^4~7| have derived the 2PN accurate gravitational waveform and the associated energy 



and angular momentum fluxes for inspiralling compact binaries through 2PN order by two 
independent methods: the BDI approach based on a mixed multipolar post-Minkowskian and 
post-Newtonian framework together with asymptotic matching and analytic continuation 
f4"8| and the recently improved Epstein- Wagoner (EW) formalism by Will and Wiseman 



f46| which provides a method to carefully handle the divergences of the older EW treatment. 
In view of the above discussion it is natural to investigate the possibility of extending 
the treatment of Iyer and Will to 2PN accuracy beyond the Newtonian (2.5PN) radiation 
reaction and this is what we propose to take up in this paper. The knowledge of the reactive 
acceleration beyond the lowest order could also have practical uses. For instance, Lincoln and 
Will |50| have studied the late-time orbital evolution of compact binaries with arbitrary mass 
ratios. They described the orbit using the osculating orbital elements of celestial mechanics 
and used the D amour- Deruelle two-body equations of motion including Newtonian radiation 
reaction terms ]2"7| , [TE| | to evolve these orbital elements. The extension of this work to include 
1PN radiation reaction is still not available. Recently, a 2PN accurate description for the 
motion of spinning compact binaries of arbitrary mass ratio was obtained in a generalized 
quasi-Keplerian parameterization initially suggested by Damour and Schafer [|5T|-|54"| . These 
orbital elements have also not been evolved to 2PN radiation reaction order. Our present 
computation is a step in that direction. These attempts to study the evolution of binary 
orbits would be complementary to those using the test particle limit f40| , |4l"| . 

To summarise: Starting from 2PN accurate energy and angular momentum fluxes for 
compact binaries of arbitrary mass ratio moving in quasi-elliptical orbits , we obtain 

the 4.5PN reactive terms in the equations of motion by an extension of the IW method. 
Schematically, the equations of motion for spinless bodies of arbitrary mass ratio are 
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a = 



^X 



mx . 



[1 + 0(e) + 0(e 2 ) + 0(e 2 - 5 ) + 0(e 3 ) + 0(e 3 - 5 ) + 0(e 4 ) + 0(e 4 ' 5 ) + 



3.5^ 



(1-1) 



where x and r = |x| denote the separation vector and distance between the bodies, and 
m = mi + 1^2 denotes the total mass. The quantity e is a small expansion parameter that 
satisfies e ~ {v/c) 2 ~ Gm/(rc 2 ), where v and r are the orbital velocity and separation of 
the binary system. The symbols 0(e) and 0(e 2 ) represent post-Newtonian (PN), post-post- 
Newtonian (2PN) corrections and so on. Gravitational radiation reaction first appears at 
0(e 2,5 ) beyond Newtonian gravitation, or at 2.5PN order. We call this the "Newtonian" 
radiation reaction. "Post-Newtonian" radiation reaction terms, at 0(e 3,5 ), were obtained by 
Iyer and Will [p8LP£ 



and Blanchet [31,32 



Here we obtain the 2PN radiation reaction, at 
0(e 4 ' 5 ). The 4.5PN reactive terms are determined in terms of twelve arbitrary parameters, 
which along the lines of |2^j29|] , are associated with the possible residual 'gauge' choice at the 
4.5PN order. These results valid for general orbits are specialized to the two complementary 
cases of circular orbits and radial infall. The expressions for r and u> for the quasi-circular 
orbits and z for radial infall to 4.5PN order are in agreement with [£3L[55| as required. We 
next examine critically the origin of the 'redundant' equations in the formalism and examine 
our understanding of this redundancy by exploring variant schemes which differ from the 
original IW scheme in their choice of the functional forms for the arbitrary terms in energy 
and angular momentum. 

The paper is organized as follows. In section 0], we describe the IW method to obtain the 



2PN reactive terms. Section [TTl] examines the question of redundant equations and explores 
'variants' of the original IW scheme that differ in their choice of the ambiguities in energy and 



angular momentum. Section [TV| discusses the question of the undetermined parameters and 
arbitrariness in the choice of the gauge, in particular at 4.5PN order. Section [V] is devoted 
to the particular cases of quasi-circular orbits and head-on infall. Section [VT| contains some 
concluding remarks. In the appendix, for mathematical completeness, we prove that the 
far- zone flux formulae and the balance equations admit more general solutions if one relaxes 
the requirement that the reactive acceleration be a power series in the individual masses of 
the binary or, equivalently, that it be nonlinear in the total mass. 



II. IW METHOD FOR REACTIVE TERMS IN THE EQUATIONS OF MOTION 

A. The Procedure 

We consider only two-body systems containing objects that are sufficiently small that 
finite-size effects, such as spin-orbit, spin-spin, or tidal interactions can be ignored. The 
dynamics of such systems is well studied and the two-body equations of motion conveniently 
cast into a relative one-body equation of motion is given by : 

a _ a i a (l) , a ( 2 ) , a (2-5) , (3) , (3.5) (4) (4) (4.5) Q( b s (l? , , 

«* — <*-N "I" dp N -f d 2 p N -f a RR -|- d 3 p N "I" <A 1RR -f a 4 p N -f d Tail -f -|-d- 2RR "1" V\t ) , I^Z.IJ 

where the subscripts denote the nature of the term, post-Newtonian (PN), post-post- 
Newtonian (2PN), Newtonian radiation reaction (RR), post-Newtonian radiation reaction 
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(1RR), 2PN radiation reaction (2RR), tail radiation reaction and so on; and the superscripts 
denote the order in e. For our purpose we need to know explicitly the acceleration terms 
through 2PN order and they are given by P7 , 55 , ^0] (G = c = 1) 
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where fi = m 1 m 2 /m is the reduced mass, with i] = fi/m, and n = x/r. The 77.5PN reactive 
accelerations are determined by following the 'What else can it be ?' procedure employed 
in IW which we summarise here. One writes down a general form for the Newtonian (e 25 ), 
1PN (e 3 5 ) and 2PN (e 4 5 ) radiation-reaction terms in the equations of motion for two bodies, 
ignoring tidal and spin effects. For the relative acceleration a = ai — a 2 , one assumes the 
provisional form 



--r,(m/r 2 )(m/r)[-(A 2 , b + A 3 , 5 + A 4 , 5 )rn + (B X5 + B 3 , 5 





B 



4.5; 



(2.3) 



The form of Eq. ( [2.3|) is dictated by the fact that it must be a correction to the Newto- 
nian acceleration (i.e., be proportional to m/r 2 ), must vanish in the test body limit when 
gravitational radiation vanishes (i.e., be proportional to 77), must be dissipative, or odd in 
velocities (i.e., contain the factors r, n and v linearly) and finally, must be related to the 
emission of gravitational radiation or be nonlinear in Newton's constant G (i.e., contain 
another factor m/r). The last condition may be more precisely stated by requiring that 
the reactive acceleration be a power series in the individual masses mi and 7772 |57| . For 
spinless, structureless bodies, the acceleration must lie in the orbital plane (i.e., depend 
only on the vectors n and v). The prefactor 8/5 is chosen for convenience. To make the 
leading term of 0(e 2 5 ) beyond Newtonian order, A2.5 and B2.5 must be of 0(e). For this 
structureless two-body system the only variables in the problem of this order are v 2 , m/r, 
and r 2 . Thus A2.5 and S2.5 can each be a linear combination of these three terms; to those 
terms we assign six "Newtonian radiation reaction" parameters. Proceeding similarly, A 3 5 
and -B3.5 must be of 0(e 2 ), hence must each be a linear combination of the six terms t> 4 , 



v 2 f 2 , 



r ;:2 



m/r, f , and (m/r) . To these we assign 12 "1PN RR" parameters. And 



4^.2 



v r 



v m/r, 

finally, A^_ 5 and -B4.5 must be of 0(e 3 ), each a linear combination of the 10 terms t> 6 
v 4 m/r, v 2 r A , v 2 (m/r) 2 , v 2 r 2 (m/r), r 6 , r A (m/r), r 2 (mjr) 2 and (m/r) 3 to which we assign 
20 "2PN RR" parameters. The 6 Newtonian RR and 12 post-Newtonian RR parameters 
were first determined in IW p8| , p9| . This solution has been checked and reproduced in the 



preliminary part of this investigation and constitutes an input to supplement the conserva- 
tive acceleration terms in Eq. ( |2.3|) for the present study. Our aim is to evaluate these 20 
parameters appearing in A 4 5 and -B4.5 that will determine the 2PN radiation reaction. It 
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is worth pointing out that in the calculation we are setting up, the terms in the equations 
of motion of 0(e 3 ) and 0(e 4 ) beyond Newtonian order do not play any role. The former is 
non-dissipative but not yet computed; the latter on the other hand includes dissipative parts 
due to the 'tail' effects |58|-|6l|] which have been separately balanced by the tail luminosity 



in the works of Blanchet and Damour |58| , |32| . However all the radiation reaction results will 
remain as 'partial results' in the saga of equations of motion until a complete treatment a 

and Damour is available through 3PN order and later through 



la Chandrasekhar 
4PN order. 

Through 2PN order, the equations of motion can be derived from a generalized La- 
grangian that depends not only on positions and velocities but also on accelerations. To 
this order, that is in the absence of radiation reaction, the Lagrangian leads to a conserved 
energy and angular momentum given by p7| , |56| , |62 | 



E — En + -Epn + -E2PN 

J = Jn + JpN + J2PN , 
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(2.5c) 
(2.5d) 
(2.5e) 

(2.5f) 



and where L^r = /xx x v. 

Through 2PN order, the orbital energy and angular momentum per unit reduced mass, 
E = E//I— |f 2 — m/r + 0(e 2 ) + 0(e 3 ), J = x x v[l + 0(e) + 0(e 2 )], are constant, and cor- 
respond to asymptotically measured quantities. However, the radiation reaction terms lead 
to non- vanishing expressions for dE/dt and dJ/dt containing the 20 undetermined parame- 
ters. Following IW, starting from the 2PN-conserved expressions for E and J we calculate 
dE/dt and dJ/dt using the 2PN two-body equations of motion p7|p6i^0| supplemented by 
the radiation-reaction terms of Eq. ( |2.3| ). In the balance approach, this time variation of the 
'conserved' quantities is equated to the negative of the flux of energy and angular momen- 
tum carried by the gravitational waves to the far-zone. Thus in addition to the EOM and 
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conserved quantities we need the 2PN accurate expressions for the far-zone fluxes of energy 
and angular momentum for a system of two particles moving on general quasi-elliptic orbits. 
The waveform, energy and angular momentum flux have been computed by Gopakumar and 
Iyer j|7| using the BDI [ 48|,45| formalism, and independently the waveform and energy flux 
by Will and Wiseman |^6[ using their new improved version of the EW |^9[ formalism. We 
quote below the final results for the fluxes per unit reduced mass: 
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In the above expressions Ljy = Lw / ' \i and the tail terms are not listed. It is important 
to emphasize that the 'tail' contribution to the reaction force is such that the balance 
equation for energy is verified for the tail luminosity f5"8",52||. This corresponds to the 'tail' 



acceleration at 4PN. With this part independently accounted for, in our analysis we focus on 
the 'instantaneous' terms without loss of generality. It is worth recalling that the 'balance' 
one sets up in the above treatment is always modulo total time derivatives of the variables 
involved. This is crucial to realize and in IW this was systematically accounted for by 
noting that at orders of approximation beyond those at which they are strictly conserved 
(and thus well defined), E and J are ambiguous upto such terms. Consequently, we have the 
freedom to add to E and J arbitrary terms of order e 2,5 , e 3,5 , and e 4,5 beyond the Newtonian 
expressions without affecting their conservation at 2PN order. There are three such terms of 
the appropriate general form at 0(e 2,5 ) in each of E and J, respectively, 6 each at 0(e 3,5 ) and 
10 each at 0(e 4,5 ), resulting in 6 additional Newtonian RR parameters, 12 additional 1PN 
RR parameters and 20 additional 2PN RR parameters, respectively. As discussed in detail 
in the following section, these numbers are very much tied up with the 'functional form' we 
assume for the ambiguous terms and in this section we follow IW in close detail. Equating 
time derivatives of the resulting generalized energy and angular momentum expressions E* 
and J* ( rather than only the conserved expressions) to the negative of the far-zone flux 
formulae and comparing them term by term one seeks to determine the extent to which one 
can deduce the 4.5PN reactive acceleration terms by the refined balance approach. 



B. The 2PN RR computation and results 

The above procedure is implemented order by order. All the computations were done 
with MAPLE |63J and independently checked by MATHEMATICA |>§. At the leading 



order, when the flux is given by the quadrupole equation, one deduces the 'Newtonian RR' 
or 2.5PN term in the acceleration. In this case, in addition to the 6 unknowns in the 
reactive acceleration, one has 3 unknowns each for the possible 2.5PN ambiguities in the 
E* and J*. As demonstrated in IW, the balance equations yield 12 constraints on these 12 
Newtonian RR parameters. Of the 12 constraints, only 10 are linearly independent, and thus 
finally one obtains 10 linear inhomogeneous equations for 12 Newtonian radiation reaction 
variables. Solving these equations one obtains explicit forms for A2.5, B2.5 and £"2.5, J2.5 in 
terms of two 2.5PN arbitrary parameters. To get the 3.5PN reactive terms, one adopts the 
above solution and extends the calculation to 0(e 35 ) after introducing £3.5 and J3.5 with 12 
additional 1PN RR parameters. At 3.5PN there are 20 constraints on the 24 post-Newtonian 
radiation reaction parameters; of the 20 only 18 are linearly independent; the solution to 
this system yields explicit forms for A 3 5 , B 3 5 and £3.5, J3.5 in terms of six 3.5PN arbitrary 
parameters. Since we need these results for the present computation, we reproduce them 
from IW |65i. 
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r r \ r / 

^ 4 9 m 9.9 .9 m .4 f m\ 2 _ „ , 

r r V r / 



where 



/1 = 2^(H7 + 13277) - |a 3 (l - 377) + 36 - 3p 5 , (2.9a 
h = -^( 297 " 310 ^) " 3/52(1 - 47/) - |a 3 (7 + 13^) - 2^ - 36 + 36 + 3p 5 , (2.9b 
/a = ^(19 - 7277) + ^3(1 - 377) - 56 + 56 + 5p 5 , (2.9c 

U = -^( 687 " 368r ?) " 6 ^ + ^«3(54 + 1777) - 26 - 56 - 66 , (2.9d 

/s = -76, (2.9e 

h = "^( 1533 + m v) - ^(14 + 977) + 3a 3 (7 + 477) - 26 - 36 , (2.9f 

gi = _ 3 (i _ 3?7 ) _ ^ 2 (i _ 377) - 6 , (2.9g 

92 = -^(139 + 76877) - |/3 2 (5 + 17 V ) + 6 - 6 , (2.9h 

53 = ^(369 - 62477) + ^(3(3 2 + 2a 3 )(l - 3 V ) + 36 - 3p 5 , (2.9i 

9a = ^(295 - 33577) + i/? 2 (38 - H77) - 3a 3 (l - H + 26 + 46 + 3p 5 , (2.9j 

95 = 4( 19 - 72r l) ~ 5«3(1 - 37/) + 5p 5 , (2.9k 
zo 

56 = "^(634 - 6677) + (3 2 (7 + 377) + 6 • (2.91 

The quantities a 3 , f3 2 , 6, 6j £3; £5 an d P5 are parameters that represent the unconstrained 
degrees of freedom that correspond to gauge transformations. In addition to the reactive 
terms listed above, one of the coefficients that determine the 2.5PN ambiguity in E and J and 
three of the coefficients that determine the corresponding 3.5PN ambiguity are nonvanishing. 
We list these also since they are needed for setting up the 4.5PN computation: 



£*1 = - (2 + p 2 ) , 


(2.10a) 


& = ~(1-4ij) , 


(2.10b) 


Ps = 6 + ^ (307 - 54877) , 


(2.10c) 


p 6 = 6- ^(271 - 21477). 


(2.10d) 
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We now adopt the 2.5PN and 3.5PN solutions given by Eqs. gTSj) . Q an d ( PTTOl) . 
Following the IW strategy, we assume the 4.5PN terms in the equations of motion to be of 
the form 



A 4 5 = h lV 6 + h 2 v 4 r 2 + h 3 v A - + h,v 2 r A + h 5 v 2 ( -) ' 

r \r J 



+h 6 v 2 r 2 — + h 7 r 6 + h 8 r A — + h 9 r 2 [ — 
r 



, m 
r 



B 



4 5 = + k 2 v 4 r 2 + k 3 v 4 — + k 4 v 2 r 4 + k 5 v 2 (— 

r \ r 



m 



, m 



+k 6 v z r z — + k 7 r b + fc 8 r 4 — + k 9 r z [ - + k 10 



my 

r / 



(2.11a) 



:2.Hb) 



We also assume for the ambiguity in £4.5 and J4.5 the restrictions and functional forms 
adopted in IW and also require that J remain a pseudo- vector. The 'generalized' "energy" 
and "angular momentum" through 4.5PN are thus given as sums of the conserved parts Eqs. 
( |2.5|) , the 'most general' 2.5PN and 3.5PN contributions - i.e., with coefficients determined 
by the Newtonian RR and 1PN RR calculations, and arbitrary 4.5PN terms. We use E* 
and J* to distinguish these quantities from the conserved energy and angular momentum. 
We get (per unit reduced mass) 



E* = E N + E PN + -E2PN + -E2.5 + -E3.5 + E. 



>4.5 



8 / m\ 2 

E N + E PN + E 2PN + -77 f —J r[(2 + (3 2 )v 2 - a 3 r 2 } 





2 






) 


r 


V r / 






/m N 


2 






) 


r 


V r / 







£it> 4 + + + e 4 r 4 + £ 5 r 2 - - ±(1 

r r 21 



ipiv + ip 2 v 4 f 2 + h ip4V jL f' 1 + 

r 



1 m 



.2-4 



m 
V r 

2-2 



/// 



+ ipQV r h 

r 



TTi \ { Tfl 

+^ 8 r q — + iJ 9 r z ( —J + V10 ( — 



(2.12a) 



1 at -r jpn -r ^2PN + J2.5 + J3.5 + J4.5 



- m 



Jn + Jpn + J2PN + -r]L N —r fo- 



rm, 



o ~ m 
--77L7V— r 

5 r 

8 ~ m 
--^Lat— t 

5 r 



1 



,m 



-(307 - 548r/ + 84&)i; a - + p 5 r 2 - - -(271 - - 42£ 3 ) - 
84 r r 42 \r 



2i 



.4-2 



.2-4 



Xiv +X2V r + X3V — + Xivr + X5V y 



m 



.2-2 



m 



+ XqV r h xir 

r 



1 m 



/ m 



+Xsf — + Xvr [ — ) + Xw ( — 



m 



(2.12b) 



We now compute the 4.5PN terms in dE*/dt and dJ*/dt using the identities 

ldv 2 



2 dt 
d(x x v) 

dt 



= v ■ a . 



= x x a . 



(2.13a) 
(2.13b) 
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where a is given by Eqs. Q, (gl|, (p|), (|pD and ( PI . To compute E* and 

J* to 0(e 4 ' 5 ), one needs to evaluate (E N ,J N ), (E 1PN , Jipn) and (E 2 pn, J2Pn) by using 
a to 0(e 4 ' 5 ), 0(e 3 ' 5 ) and 0(e 2 ' 5 ), respectively. On the other hand, for time derivatives of 
the 'ambiguity parts', (-E4.5, J4.5), (-E3.5, J3.5) and (£"2.5, J2.5), the relevant accelerations 
are the 'conservative' accelerations to order Newtonian, post-Newtonian and second post- 
Newtonian, respectively. Schematically, we get, 



dE* 
~~dt 



dJ* 
~dT 



m 
15 r 6 



2 r 



m 
r 



12v 2 
m 



llr' 



m 
r 



1 

28 



2-2 



(785 - 852t])v 4 + 2(-1487 + 1392rj)v 2 r 
m 



+ 160(-17 + 77) — v 2 + 3(687 - 620ry)r 4 + 8(367 - 15^)— r 



+ 16(1-4^) ( 



m 
V r 



15 



i=l 



8 ~ m 
5 r z 



m 
r 



2v l + 2- 



m 



3r 2 ) + 



m 



1 

84 



(307 - 548t/)?/ + 6 (-74 



-4(58 + 9577)— v 2 + 3(95 - 360r?)r 4 + 2(372 + 197t/) — r 2 



+2 (-745 + 277) 
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+ E< ? P 1 3? 1 



i=l 



(2.14a) 
277^) uV 

(2.14b) 



where 



1 ... 15) 



r / 



4.4 4 

v r ,v 



m 

i" I — 

r 



r / V r J \ r . 





•2 2 / m 


(7) 


r ,iT — 




V r 






) ^ 


(")> 







7A 6 , 



■,6 



m 

— I r", r" 
r 



(2.15) 



and 7?] 4 ' 5 ' and iSj 4 ' 5 ^ consist of combinations of the parameters hi and fcj from A 4 5 and 
-^4.5; ^i) Xi combined with functions of 77 from E i5 and J4.5, £1, £2, £3, £4, £5, Ps combined 
with functions of 77 from 1PN corrections of 3.5PN terms and 03 and 02 combined with 
functions of rj from 2PN corrections of 2.5PN terms. We equate dE* /dt and dJ*/dt thus 
obtained to the negative of the 2PN far-zone fluxes given by Eqs. (|2.7| ). This results in 
30 constraints on the 40 parameters hi, ki, ipi and \%- Two of these constraints being 
redundant, of the 30 constraints only 28 are linearly independent. The system of 28 linear 
inhomogeneous equations for 40 variables is therefore under-determined to the extent of 12 
arbitrary parameters, and we choose these to be ipi . . . ipg, Xe, Xs an d X9- With this choice, 
the coefficients in Eq. ( |2.11| determining the 4.5PN reactive acceleration are given by 



hi 



168 

+3^ 2 



121 - 22787/ + 40 127/' 
3xe , 



Wl - 9t? + 2177 2 ) - -(£ 2 - p 5 )(l - 377) 



(2.16a) 



^(329 - 1487r/ + 12447/ 2 ) + ~a 3 (l - 9t/ + 21 V 2 ) + ^ 2 - & - p 5 )(l - 3 V ) 
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-5^2 + 5^4 + 5xe - 5xs , (2. 16b) 

h 3 = ^(7692 - 874297/ + II2I877 2 ) + -a 3 (l - 97r] + 25t7 2 ) + ^/3 2 (3 - 3r] - 19r] 2 ) 

+36(1 - 477) - |(& - p 5 )(7 + 1377) - ^5(1 - 377) 

-2^i - 3^ 2 + 3^ 6 + 3xe - 3X9 , (2.16c) 
= -^(39 - 16377 + 97T7 2 ) + ^ 4 (1 - 377) - 7^ 4 + 7^ 7 + 7x 8 , (2.16d) 

h 5 = --L(37089 - 640057/ + 1129777 2 ) + 9a 3 (2 + 13r? + 2r] 2 ) + i/? 2 (48 - I2I77 - 54t7 2 ) 
252 4 

+6(14 + 977) + 3(6 - p 5 )(7 + 477) + 36(1 - 477) - ^6(7 + 1377) 

-2^3 - 3t/; 6 + 3t/; 9 + 3x 9 , (2.16e) 

h 6 = _JL(45475 - 21953577 + 4312177 2 ) - 7« 3 (14 - 40377 + 77t7 2 ) - -#277(7 - 137/) 
504 4 2 

+6r/6 + ^2(68 - 977) - ^ 4 (7 + 1377) + 36(1 - 377) - ip 5 (62 + I977) 

-4^2 - 5^4 - 6t/> 6 + 5t/; 8 + 2 X6 + 5x 8 + 6x9 , (2.16f) 

hT=-tol> 7 , (2.16g) 

1 1 33 

/i 8 = ^( 5002 - 36589r ? + 4496T7 2 ) - -a 3 7?(233 - 63t/) + — p 2 r)(l - 3r)) + 3t/6 
252 8 4 

+ ^6(82 + 2377) + 577p 5 - 27>4 - 7^ 7 - 8^ 8 , (2.16h) 

h 9 = _L(181371 - 3424797/ + 42598T7 2 ) - -a 3 (117 + 1097/ + 6t7 2 ) 

-^ 2 (28 + 24577 + 2077 2 ) + 2776 + (26 + 56)(7 + 477) + 7776 + ^6(60 + 2177) 
+377P5 - 2^ 6 - 57> 8 - 7t/> 9 , (2.16i) 
h w = ^(265265 + 2622307/ + 1507277 2 ) - -a 3 (102 + 1777/ + I677 2 ) 

+ ^ 2 (200 + 32577 + 4O77 2 ) + 6(14 + 977) + 36(7 + 477) - 2t/> 5 - 3t/> 9 , (2.16j) 

h = + 2)(1 - 77 - II77 2 ) + ^6(1 - 37/) - , (2.16k) 

k * = -T7^( 499 " 265677 - 1467/ 2 ) - ^a 3 (l - 3rj - 3t7 2 ) - |/3 2 (1 - 77 - II77 2 ) 
loo 2 o 

-^(36 - 26 - P5)(l - 377) + 3^ - 3X6 , (2.161) 

k z = 7^7(81 - 912777 - 1448277 2 ) - |/3 2 (3 + I2I77 + 7 V 2 ) + ^(5 + 17tj) 
oU4 o 2 

+ 2&(1 - 377)+7/> 1 -t/>3, (2.16m) 

h = ^(329 - 148777 + 124477 2 ) + |a 3 (l - 377 - 3 V 2 ) - ^(2& - 26 + p 5 )(l - H 

+5X6-5X8, (2.16n) 
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h = -r^r(H07 - 80577 - 50877 2 ) + j/3 2 (16 + 255t/ + 22t/ 2 ) - 6(7 + 377) 

+^3(5 + 1777)+^-^, (2.16o) 

fc 6 = 77Tt(1797 + 548167/ - 22463t/ 2 ) + -a 3 (l + 3ry + 5t/ 2 ) - ^/3 2 (42 - 485t/ + 173t/ 2 ) 
5U4 Z 4 

_ 1^(56 _ 4977) - 3(6 + 26 - 6)(1 - 3t/) + ^p 5 (7 + H77) 

+4^ + 4^ + 3x6-3x9, (2.16p) 
k 7 = -^(39 - 1637/ + 97t/ 2 ) - 76(1 - 3t/) + 7 Xs , (2.16q) 

A; 8 = — ^-(39808 - 92788t/ + 24563t/ 2 ) + -a 3 (U - 1057/ + 59t/ 2 ) - -/? 2 t/(69 + 13t/) 
504 2 8 

-3t/6 - (26 + 56 + 66)(1 - 3t/) - ip 5 (62 + 3t/) + 2 Xe + 5 X s + 6 X 9 , (2.16r) 

A; 9 = —(83 1 9 - 76837/ + 11809t/ 2 ) + 3a 3 (3 - 13t/ - rf) - -/? 2 (194 + 215t/ + 24t/ 2 ) 
252 4 

-(26 + 3p 5 )(7 + 3t/) - ^6(44 - 9t/) - 36(1 - 3t/) + 2t/> 3 + 5^ 5 + 3 X9 , (2.16s) 
fc io = 7^^(425413 + 1116367/ - 6912t/ 2 ) - J/3 2 (53 + 103t/ + 4t/ 2 ) 

ZZOo Z 

-6(7 + 3t/)+t/> 5 . (2.16t) 

At the 4.5PN order, 4 parameters determining £4 5 and J4 5 are non-vanishing and are given 
by 



^10 = 7^77 (362 - 15487/ + 400t/ 2 ) , 



189 

Xs = fa + 7777 (2665 - 123557/ + 12894t/ 2 ) , 
504 v ' 

X5 = ^3 + ^ - 77^77 (524 - 44837/ + 3675t/ 2 ) , 



X10 = ^5 - ^ + 7^ (775 - 39397/ + 2942t/ 2 ) . (2.17) 

A final minor remark is with regard to the two possible ways one may implement the require- 
ment that the ambiguity in J* be a pseudovector. One may either choose it proportional to 
Lin as in the treatment above or to the conserved angular momentum J. At 2.5PN order 
both choices are identical. At the 3.5PN order, the two choices lead to an identical system 
of linear equations barring a translation in the values of p 3 and p 6 by an amount given by 
the coefficients of v 2 and m/r in Jipn- 

Ps^Ps = Ps + 7j(l — 3?7) /3 2 , 

P6^P"6 = P6 + (3 + 7/) fo . (2.18) 

Since p 3 and p$ are not among the arbitrary parameters determining the solution, the solution 
determining the reactive terms and 6 is unchanged. Only the expressions for p 3 and p 6 are 
changed to 
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f>3 = 6 + ^ (307 - 548?/) + ~ (1 - 377) (3 2 , 

Pe = 6 - ^(271 - 21477) + (3 + 77) ft ■ (2.19) 

At 4.5PN order, however, the situation is different. Indeed, as before, the two choices lead 
to an identical system of linear equations barring a translation in the values of the five 
parameters X3, X5, Xe, X9 and Xw- 

Xs^Xs = X3 + l(l~Qv + 21r/ 2 ) & - i (1 - 377) 6 - ^ (307 - 14697? + 1644r/ 2 ) , 
X5->Xis = Xs + J (l + 677 - 3t/ 2 ) ft - (3 + 77) £1 - ~ (1 - 377) £3 



2 v , , 2 
-—(325- 15577 - 595T7 2 

Xe^Xe = Xe - g ( 1_ 3r/) p 5 , 
X9^X9 = X9 - ^ (2 + 577) 77ft - (3 + 77) p 5 , 
Xio-Xio = X10 - \ (22 + 65t/) ft - (3 + 77) £ 3 + ^ ( 5691 " 2597 ^ " 1498 ^ 2 ) ■ ( 2 - 20 ) 

Consequently, in terms of the above 'shifted' variables, the solutions for the reactive accel- 
erations are identical. As xe an d X9 are among the independent parameters that determine 
the reactive acceleration, in terms of x& an d X9 the two choices yield equivalent but different 
looking solutions for the 4.5PN reactive terms in the equations of motion. 

Of the two choices, the second choice is more convenient for calculations by hand since 
dJ/dt = to 0(e 2 ), but has no special advantage when the calculation is done on a computer. 



III. REDUNDANT EQUATIONS AND RELATED VARIANT SCHEMES 

It was noticed in IW that both at the 2.5PN and at the 3.5PN order, the 'balance 



procedure' leads to two redundant constraint equations p9fl . Here, at 4.5PN order, we once 
again obtain two redundant constraint equations. In this section, we examine critically the 
origin of these redundant equations. 



In implementing the 'refined balance procedure' for the general orbits, IW |U balance 
the 'energy flux' and 'angular momentum flux' completely independently of each other. 
However, for circular orbits, these fluxes are not independent but related |66| via: 



d£\ 
~di) 



v 2 j 



far— zone 



where J is defined by the equation 

/WT\ 



~dt) 



far— zone 
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The general balance should reflect this limit and we find that for Newtonian RR a linear 
combination of the 6 equations representing energy balance and another linear combination 
of the 6 equations representing angular momentum balance are indeed identical and given 
by: 

ei + e 2 -4 = 0. (3.1) 

Similarly at 3.5PN we have 

9i + 92 + 96 ~ (3 - rj) (3 2 + ^ (2927 - 252 V ) = , (3.2) 
and finally at 4.5PN order the 'degenerate' equation is 

h + k 3 + k 5 + kw + ( 3 ~ ^ + ^ + ^( 90 + 13r ? + 6r ? 2 )& 

— (635771 + 297117r/ - 8IOOO77 2 ) = . (3.3) 

4536 

Thus we can trace the existence of one of the redundant equations in the IW procedure to the 
fact that for circular orbits the energy and angular momentum fluxes are not independent 
but proportional to each other. 

The mystery of the other redundant equation was not so easy to resolve but after a careful 
examination of the system of equations and 'experiments' in modifying the system, we could 
finally track it back to its source. The observation that this redundant equation relates 
the coefficients of the polynomial representing the ambiguity in J led us to examine the 
functional form that IW proposed as the starting ansatz for the calculation. A comparison 
of the functional forms for the ambiguity in E and J Eqs. ( (|2.12| )) reveals that indeed IW 
assume a more general possibility for J than required. The ambiguity in angular momentum 
leads to terms more general than required by the far-zone flux formula and time derivative of 
the leading term using the reactive acceleration. The absence of such terms in the far-zone 
flux then yields only the trivial solution for these additional variables in J, and the second 
redundant equation is just a homogeneous linear combination of these trivial solutions. Thus 
the second redundant equation in the IW scheme is due to the fact that the IW scheme - 
extended here to 4.5PN order — is not a 'minimal' one. 

To verify this 'conjecture' we experimented with alternatives for the functional form that 
one assumes as the starting expression for the ambiguity in E and J — the 2.5PN, 3.5PN 
and 4.5PN order terms. In the first instance, we replace the IW scheme — labelled for 
clarity of reference by IW21 — by the 'minimal' variant in Eq. ( |2.12j ) — labelled by IW22. 
The notation IW21 indicates e.g., that (m/r) 2 is pulled out in E while only (m/r) 1 is pulled 
out in J. As explained above, the minimal choice for J* is obtained by pulling out the factor 
(8/b)rjh^{m/r) 2 r from arbitrary terms in J*, rather than the factor (8/b)rjh^{m/r)r as in 
the IW scheme for J*. This reduces by one the order of the polynomial in v 2 , r 2 , and m/r 
that constitutes the arbitrariness, and consequently implies a reduction in the number of 
variables that characterise the ambiguity in J to one for J 2 .s, three in J 3 5 and six in J 4 5 . 
Thus in the IW22 scheme, at the 2.5PN level we have 6 variables in the reactive acceleration, 
3 variables determining the energy ambiguity E2.5 and 1 variable determining the ambiguity 
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in J 2.5 i.e., 10 variables in all. The balance equations lead to 9 equations — 6 from energy 
and 3 from angular momentum — of which 8 are linearly independent. In other words, 
there is only one redundant equation. The linear system of 8 equations for 10 variables is 
then the same as before and leads to the IW21 solution in terms of 2 arbitrary parameters. 
(The two extra variables in IW21 are identically zero.) Similarly, at the 3.5PN level we 
have 12 variables in the reactive acceleration, 6 variables determining the energy ambiguity 
E 3 5 and 3 variables determining the ambiguity in J 3 5 , i.e., 21 variables in all. The balance 
equations lead to 16 equations - - 10 from energy and 6 from angular momentum — of 
which 15 are linearly independent, leaving only one redundant equation. The linear system 
of 15 equations for 21 variables is then the same as before and leads to the IW21 solution 
in terms of 6 arbitrary parameters. (The three extra variables in IW21 are identically zero.) 
Finally, at the 4.5PN level, we have 20 variables in the reactive acceleration, 10 variables 
determining the energy ambiguity -E4.5 and 6 variables determining the ambiguity in J 4 5 , 
i.e., 36 variables in all. The balance equations lead to 25 — 15 from energy and 10 from 
angular momentum — equations of which 24 are linearly independent, again leaving only 
one redundant equation. The linear system of 24 equations for 36 variables is the same as 
before and leads to the solution obtained in the previous section in terms of 12 arbitrary 
parameters. (The four extra variables in the IW21 scheme are identically zero.) The IW22 
(minimal) scheme thus confirms the conjecture that the occurence of the second redundant 
equation is special to the IW scheme (IW21) and is related to the choice they make for 
the functional form of the J ambiguity by pulling out only one factor of nonlinearity m/r 
rather than its square — the minimal choice. To double check the above explanation, we 
performed another experiment by examining a variant that would generate an increased 
number of redundant or degenerate equations. This scheme denoted by IW11 differs from 
IW21 in that the ambiguity in E* is assumed to have (8/5)r](m/r)r as the common factor, 
i.e., by pulling out only one order of nonlinearity m/r rather than its square as in IW21; 
the polynomial representing the ambiguity in E is consequently of one order more than in 
IW21. In this case, at 2.5PN order one has 6 + 6 + 3 = 15 variables and 10 + 6 = 16 equations 
of which 3 are redundant. The 13 equations for 15 variables thus yield the required solution 
in terms of 2 arbitrary parameters and similarly for higher orders. One may also explore 
the most general of choices in which only (8/5)77 ls pulled outside and the ambiguity is the 
highest order polynomial consistent with the order of the approximation. We studied one 
such scheme (IW00) in the Newtonian RR case. For convenience, the various experiments 
are summarized in Table |I[ 

To conclude: at 2.5PN, 3.5PN and 4.5PN orders all variants of IW examined in this 
subsection with different forms of the ambiguities in E and J — minimal (IW22) or IW11 
- lead to identical reactive accelerations including their gauge arbitrariness. 

At this juncture one may wonder about the issues of the 'uniqueness' and 'ambiguities' 
of the schemes discussed earlier. In this regard, we would like to make the following general 
remarks. For general orbits, in addition to the balance of energy one must take into account 
the balance of angular momentum. Thus, schemes involving only energy balance are not 
relevant except in special cases like 'circular orbits' and 'radial infall' (see section 0). Can one 
have schemes where one implements both energy and angular momentum balance but does 
not take into account the possible ambiguities in E and J? One can show that even at the 
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2.5PN level this system of equations is inconsistent! Further, is the ambiguity necessary both 
in E and J ? If one examines a scheme with both energy and angular momentum balance 
taking account of the ambiguity only in E one does obtain a consistent solution upto 4.5PN 
order but with only half the number of arbitrary parameters as in the IW scheme. The 
reduced 'gauge' freedom is not adequate to treat as special cases the Burke-Thorne gauge 
at the 2.5PN level or the Blanchet choice at the 3.5PN level. And finally, in a scheme with 
both energy and angular momentum balance taking account of the ambiguity only in J one 
obtains a consistent solution at 2.5PN order containing no arbitrary parameters at all. No 
solution is possible at higher orders. 

On general considerations, the reactive acceleration should be a power series in the 
individual masses mi and m 2 or equivalently, it should be nonlinear in the total mass m as 
assumed in earlier sections. It is interesting to investigate whether the functional forms of 
the far-zone fluxes and the balance procedure necessarily lead to such 'physical' solutions 
alone or whether they are consistent with more general possibilities. In the Appendix, for 
mathematical completeness |)7] we investigate this question in detail and prove that the flux 
formulas and balance equations do not constrain the reactive acceleration to their 'physical' 
forms alone but allow for a more general form for the reactive acceleration. 



IV. ARBITRARINESS IN REACTIVE TERMS AND GAUGE CHOICE 

It is well known that the formulas for the energy and angular momentum fluxes in the 
far- zone are gauge invariant, i.e., independent of the changes in the coordinate system that 
leave the spacetime asymptotically flat. On the other hand, the expressions for the reactive 
force are 'gauge dependent' and consequently e.g., the Chandrasekhar form is different from 
the Burke-Thorne or Damour-Deruelle forms. In IW it was shown that the Burke-Thorne 
gauge corresponds to the values /3 2 — 4 and «3 = 5, while the Damour-Deruelle choice 
corresponds to /3 2 = — 1 and «3 = 0. It was further shown that the reactive acceleration 
implied by Blanchet's first principles determination of the 1PN radiation reaction indeed 
corresponds to a particular choice of the arbitrary parameters in the IW solution. One of the 
satisfactory aspects of IW was the demonstration that the part of the reactive acceleration 
not determined by the balance requirement was precisely related to the possible ambiguity 
in the choice of the gauge at that order. (The flux is equal to the time variation of the 
conserved quantities only upto total time derivatives; this ambiguity may be absorbed in a 
'change' in the relative separation vector as discussed below.) 

Following IW, we seek to establish the correspondence between the arbitrary parameters 
contained in the radiation reaction terms and the residual gauge freedom in the construction. 
The residual gauge freedom arises from the fact that the far-zone fluxes Eqs. ( |2T6| ), (|2.7|) are 
independent of changes in the coordinate system that leave the spacetime asymptotically 
flat. These coordinate changes will induce a change in x which is the difference between the 
centers of mass of the two bodies Xi(t) and x 2 (t) at coordinate time t. Following IW, we 
choose the transformation to be of the form x —>■ x' = x + <5x, where 5x can depend only on 
the two vectors x and v, 

^X = (/ 2 .5 + /3.5 + h.5)rx + (#2.5 + 93.5 + #4.5)?"V . (4.1) 
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In order that <5x/x be 0(e 2 5 ), 0(e 3 5 ) and 0(e 4 5 ), / 2 .5 and g 2 .5 must be 0(e 2 ), / 3 . 5 and g 3 . 5 
must be 0(e 3 ) and / 4 . 5 and #4.5 must be 0(e 4 ). As for the other variables, the /'s and g>'s 
will also be polynomials in the variables m/r, v 2 and f 2 . As pointed out in |[29|| , we do 
not independently take into account changes in the coordinate time t since the v-dependent 
term in 5x includes this contribution via x(t + St) ~ x(t) + vSt. 

In |29| it was proved that to cancel the dependence on the two 2.5PN arbitrary parameters 
and the six 3.5PN arbitrary parameters, Sx. should be chosen such that 
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where P a bS and Q a b's are given by 

P21 = 



2 1 

£2+ g&-P5 - 2 a s(l — 377) 



-P 2 2 
P23 
Q21 
Q22 
Q23 



6 + 6 - 2& - Ps - ^ + ^« 3 (4 + llr/) 



1 

6 + ^6 + ^6 + ^(3/3 2 - 2a 3 )(l - 37?) 



3 63 1 

66 + 56 - 36 + 56 - -6 + Ps - y/? 2 r7 - -a 3 (4 - 55/7) 

:6 + P5 -a 3 (l -3/7) 



We provisionally choose the 4.5PN part of Sx to be of the form 



h, = - 5 v 
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) 2 + P 45 (™) r 2 + P 46 r 4 



,m 



Q41W + Q42W — + V + Q 44 ( — ) + Q45 ( — )'"+ <24 6 r 



m 



■ 1 



(4.2a) 
(4.2b) 
(4.2c) 
(4.2d) 

(4.3a) 
(4.3b) 
(4.3c) 
(4.3d) 
(4.3e) 
(4.3f) 

(4.4a) 
(4.4b) 



The change in the 2PN equations of motion Eqs. ( |2.2| ) produced by this change of variable 
Eq. (44) can be determined using the known form of Sx upto 3.5PN order Eqs. ( |4.2|) , ( |4.3| ). 
the provisional form chosen above for the 4.5PN terms Eq. ( |4.4| ) and the transformations 
given below: 

x— >x' = x + <5x , 

dx dSx 
v^v =v + ^v=- + — , 
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2 12 
V — >V 



x <5x pn 



v 2 + 



2v- 



n ■ Sx) , 



eft 



rr + <5x • v + x 



d<5x 



(n ■ 5x)r 



(4.5) 



The gauge change generates reactive terms and the requirement that this change should 
cancel the dependence of the radiation-reaction terms on arbitrary parameters dictates that 
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377) 
(4.6b) 
(4.6c) 



(4.6d) 

(4.6e) 
(4.6f) 

(4.6g) 



(4.6h) 
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Q43 



1(1 -3?7- 3r/ 2 ) - ^(26 + 2U + p 5 )(l - 3^) 




2 , 16 , 1 2 
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(4.61) 



The above computation shows that as at the 3.5PN order the (12 parameter) arbitrariness 
in the 4.5PN radiation reaction formulas reflects the residual freedom that is available to 
one in the choice of a 4.5PN accurate 'gauge'. Every particular 4.5PN accurate radiation 
reaction formula should correspond to a particular choice of these 12 parameters. 

V. PARTICULAR CASES: QUASI-CIRCULAR ORBITS AND HEAD-ON INFALL 

In this section we specialise our solutions valid for general orbits to the particular case 
of quasi-circular orbits and radial infall and verify that they indeed reproduce the simpler 
reactive solutions one would obtain if one formulated the problem ab initio appropriate 
to these two special cases. We first consider the quasi- circular limit that is of immediate 
relevance to sources for the ground based interferometric gravitational wave detectors. In 
this particular case, the reactive acceleration may be deduced using only the energy balance. 
Using the reactive acceleration we compute the 4.5PN contribution to r and Co. We also 
discuss the complementary case of the radial infall of two compact objects of arbitrary mass 
ratio and determine the 4.5PN contribution to the radial infall velocity for the two special 
cases: radial infall from infinity and radial infall with finite initial separation. 



Using our general reactive solution we can compute the physically relevant quantities r 
and Co for quasi-circular inspiral, where r and 00 are the orbital separation and the orbital 
angular frequency in harmonic coordinates, respectively. As would be expected, these results 
are independent of the arbitrary parameters that are present in the reactive solution. We 
obtain the radiation reaction contribution to a upto 4.5PN for quasi-circular inspiral by 
setting r = + 0(e 2 - 5 ) and using 



A. Quasi-circular inspiral 
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in Eqs. (H), Q, ([D3 



) and (|2~T6|) . We get 
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(5.2) 



It is worth noting that for quasi-circular inspiral the energy flux determines the reactive 
acceleration without any gauge ambiguity. All the arbitrary terms in energy are proportional 
to r and hence play no role in this instance. Inverting Eq. (|5.1|), we get 
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1 + (3 - 7])v 2 + -(48 - 89r/ + Arf )v A 



(5.3) 



Differentiating Eq. ( |5.3| ) w.r.t t and noting that the a that appears is the total acceleration 
( conservative + reactive) we get, after some rearrangement 
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Using Eq. ( |5.4| ) and the expression for angular velocity {u = v/r) 
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we may express u as 



96 



ijj- 



r]{muj) 



2 /743 11 
[ma;) a ( -z^z + —r) 



+ 



/34103 13661 



V336 4 
59 



V 18144 



77 H 77* 

2016 ' 18 ' 



(ma;) 3 



(5.6) 



The results Eqs. (|5.4j ) and (]5.6|) are in agreement with |43] as expected and required, 
suggesting that the reactive terms obtained here could be used to evolve orbits in the more 



general case also [58 



B. Head-on infall 



Recently Simone, Poisson and Will have obtained to 2PN accuracy the gravitational 
wave energy flux produced during head-on infall and starting from these formulas one can 
deduce ab initio the reactive acceleration in this limit adapting IW to the radial infall case. 
As required, these results match exactly with expressions obtained by applying radial infall 
limits to the general orbit solutions and we summarize the relevant formulas in this limit in 
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what follows. Equations representing the head-on infall can be obtained from the general 
orbit expressions by imposing the restrictions, x = zn, v = in, r = z and v = r = z. For 
radial infall the conserved energy Eq. ( |2.5| ) to 2PN order then becomes 
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(9 + 7?7 + 8r] 2 ) rfz 2 (2 + 1577) 7 3 



(5.7) 



where 7 = m/z. Unlike the quasi-circular inspiral, for head-on infall we can distinguish 
between two different cases. Following [^] we denote them by (A) and (B), respectively, 
and list the expressions relevant for our computations. In case (A), the radial infall proceeds 
from rest at infinite initial separation, E(z) = E(oo) = 0, and inverting Eq. (|5.7Q we get 
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+ r [is 

In case (B), the radial infall proceeds from rest at finite initial separation z$, which implies 
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We obtain as in case (A), an expression for i given by 
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(5.10) 



where 70 = tu/zq. We first compute the 4.5PN contribution to i for case (B),the radial 
infall from finite initial separation. We use the radial infall restriction along with Eq. (|5.10 ) 
in Eqs. @, 0, ( pH] ) and (|27T6| ) to obtain 4.5PN terms in z as 
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To obtain the 2PN reactive terms for case (A), the radial infall from infinity, we use in 
Eqs. (|2.3|) , (|2.8|) , (|2.11|) and ( |2.16|) the radial infall restriction and Eq. (|5.8|) . The expression 
thus obtained is the same as obtained by putting 7 = in Eq. ( |5.11| ). The ( 's in Eq. ( |5.11| ) 
are given by 

Ci = a 3 - (3 2 , 

( 2 =Zl + t2 + U, 

C3 = £3 + £5 5 

C4 = ^3 + tp6 + ^8 , 

C5 = ^1 + ^2 + ^4 + V>7 , 

C 6 =V>5+V>9- (5-12) 

We have also computed the 2PN reactive terms for cases (A) and (B) ab initio using the 
IW method adapted to radial infall. In this case, only energy balance is needed as J = for 
head-on infall. The result thus obtained is in agreement with Eq. ( |5. 11 ). Eq. (|5.11|) may be 
integrated straightforwardly to obtain the 4.5PN contribution to z 2 in case (B) and it yields 



77o 



315 



7o 2 



16(2 7 -2 7o )^|l_ A 
5 l2r s; ' 105 

(^(-18054 + 1323177) + 1(438 - 331//), , - 2, 2 - „ ) 



+ 
+ 



(202 (1926 " 7597r?) + iis (660 ~ 25Mr])Cl + 2Cs )"° ~ 



("315 



^(-342 + 34177) + 1^(240 - 2807/)C 1 ) 7 2 7 



+ (9l5 ( - 684 + 682 ^ ) + 47k 



+ 



(1440- 1680r/)Ci) 7 o 
/ 1091065 5649317/ 5258809t/ 2 1 



V 712£ 



2079 



66528 



352 



(27156 - 498167/ + 15057t/ 2 )G 



+ 1(7 66 _ 5 27 7/)C2 + ^(546 - 417r/)C 3 - 2( 4 - 4( 5 - Ce j " : 



24 



, 21548237 26019487,, _ 275038V + _i_ (26316 _ + 6723 2) 

V 224532 49896 11088 528 v ' ' ; 

-^(4416 - 6241?7)C 2 - ^(546 - 2023r ? )C 3 + 2( 4 + 8(5)707" 

/3823453 _ 1^81430, 4399627t7 2 _ 1 _ ^ 

V 149688 14553 22176 2464 v ' ' ;s 

+^(53262 - 20274l77)C 2 + ^(24 - 28 V )( 3 - 4C 5 )7oV 

/ 567739 60899277 228227t7 2 1 , 2 , > 

+ + 504 + IO8O77 - 162277 2 Ci 

V187110 72765 27720 385 V ' MS 

1 , o s , 1 , x> \ <, Z567739 1 21798477 

(1056 - 1232t ? 2 )C 2 + —(96 - 112t ? )C 3 7o7 + 7^7^ + 



2541 v ' ^ 385 v iwjiui ^280665 218295 

228227T7 2 1 1 
~~ ! " (1008 + 2160r7 - 3244t] 2 )Ci - ———(6336 - 7392t7)C 2 



41580 1155 v ' ' ^ 22869 

1 (10560 - 12320t7)C 3 )7o 



63525 



(5.13) 



We obtain the 4.5PN contribution to i 2 for case (A) by putting 70 = in Eq. ( B.13| ). Unlike 
in the case of quasi-circular inspiral the expressions in the head-on or radial infall cases are 
dependent on the choice of arbitrary variables or the choice of 'gauge'. 



VI. CONCLUDING REMARKS 

Starting from the 2PN accurate energy and angular momentum fluxes for structureless 
non-spinning compact binaries of arbitrary mass ratio moving on quasi-elliptical orbits we 
deduce the 4.5PN reactive terms in the equation of motion by an application of the IW 
method. The 4.5PN reactive terms are determined in terms of twelve arbitrary parameters 
which are associated with the possible residual choice of 'gauge' at this order. These gen- 
eral results could prove useful to studies of the evolution of the orbits. The limiting and 
complementary cases of circular orbits and head-on infall have also been examined. 

We have systematically and critically explored different facets of the IW choice like the 
functional form of the reactive acceleration and provided a better understanding of the origin 
of redundant equations by studying variants obtained by modifying the functional forms of 
the ambiguities in E* and J*. The main conclusions we arrive at by this analysis are 

• In terms of the number of arbitrary parameters and the corresponding gauge trans- 
formations, the IW scheme exhibits remarkable stability for a variety of choices for 
the form of the ambiguity in energy and angular momentum. The different choices 
merely produce different numbers of degenerate equations. This indicates the essential 
validity and soundness of the scheme. These solutions are general enough to treat as 
special cases any particular solutions obtained from first principles in the future. 

• Relaxing the requirement of nonlinearity in m or more precisely the power series be- 
haviour in m! and 7772 permits mathematically more general solutions for the reactive 
accelerations involving more arbitrary parameters. Solutions more general than the 
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ones discussed in the appendix, e.g., a solution involving 6 parameters at the New- 
tonian level, cannot be gauged away either by gauge transformations of the form 
discussed by IW or by more general gauge transformations that differ in their powers 
of nonlinearity (jn/r dependence). However, none of these solutions are of 'physical' 
interest to describe the gravitational radiation reaction of two-body systems. 
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APPENDIX A: THE GENERAL SOLUTION TO THE BALANCE METHOD 

1. The 2.5PN reactive solution 

It should be noted that all the discussion in section |IT| follows only after one has assumed 
a functional form for the reactive acceleration — in particular, the intuitive requirement 
that it be nonlinear, i.e., contain an overall factor of m/r. It is pertinent to ask whether 
more general possibilities obtain, consistent with the far-zone fluxes, if one relaxes this 
requirement. We have explored this question in detail at the 2.5PN level and we summarize 
the results in what follows. In this instance the reactive acceleration is assumed to be: 

A2.5 = a[v* + a' 2 v 2 ™ + a' 3 v 2 r 2 + a' 4 (^) + a' 5 C^j r 2 + a' & r A , 

£ 2 . 5 = b\v A + b' 2 v 2 ™ + b' 3 v 2 r 2 + b> 4 (™) 2 + b' 5 (™) f 2 + 6 6 r 4 , (Al) 

i.e., it is determined by 12 reactive coefficients instead of the earlier 6. Recall that the 
nomenclature IW22, IW21 and IW11 refers to the functional forms chosen for the ambiguity 
in energy and angular momentum and we introduce similar notation EJ22, EJ21 and EJ11, 
respectively, in this appendix, where the acceleration has a more general form as given by 
Eq. ( |A1|) . With this form of the reactive acceleration, however, one gets e.g., in the EJ21 
scheme at 2.5PN 

E*= E N + E 2 .5 

~ 8 / m\ 2 . , 9 m . 9 . . , „ . 

= E N - -7} — r(a l v z + a 2 h a z r z ) , (A2a) 

5 V r J r 

J*= Ltv + J2.5 

— Pi — 777 777 

= L N + - V L N -r(/3 lV 2 + (3 2 - + (3 3 f 2 ) , (A2b) 
or r 
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The derivatives of E* and J* with the new form of the reactive acceleration are given by 
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Using Eqs. (|A2|) and ( |A3| ) one can understand the counts of the various variables summarised 
in Table 0. 

One can explain the new counts for the arbitrary parameters by comparing e.g., the EJ21 
scheme with a general form for the reactive acceleration as in this section with the IW21 
scheme with the restricted form for reactive acceleration as in section |TJ. One has 6 extra 
variables and 4 extra equations. However one gains an extra equation because one of the 
degeneracies is lifted. The resulting 5 equations for 6 variables lead to an extra arbitrary 
parameter resulting in a 3 parameter solution in this instance. All the other entries in Table 
H can be similarly understood by comparison of Tables [I] and |I|. 

The reactive solution resulting from the EJ22 scheme in this instance is exactly the same 
as the IW21 reactive solution discussed earlier. From the EJ21 scheme one obtains a solution 
with three arbitrary parameters given by 
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This construction can be generalised to 3.5PN and 4.5PN orders in which cases the number 
of arbitrary parameters are 8 and 15, respectively. The EJ11 and EJ00 schemes on the other 
hand lead to a solution with six arbitrary parameters at the 2.5PN level. However, not all 
these solutions are similar in regard to the possibility of gauging away all the arbitrary 
parameters they contain. 



2. The 2.5PN gauge arbitrariness 

We have also investigated the question whether all the extra arbitrary parameters ap- 
pearing in schemes with the general form of reactive acceleration (See Table f|.) can be 
gauged away? We find that at 2.5PN order, though this is possible with the 3 parameters 
of the EJ21 scheme, it is not true for the 6 arbitrary parameters in the EJ11 and EJ00 
schemes. For this reason the EJ11 and EJ00 schemes are not satisfactory and we discuss 
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them no further. We present here for the EJ21 scheme details of the gauge calculation at 
2.5PN order. We choose 5x to be 



5x = y (7) (As™ + #2.5 rV ) > ( A5 ) 



where 5 and g' 2 5 are given by 
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For the reactive acceleration given by Eqs. (|A1|) and (|A4|) we obtain 

(A7a) 

(A7b) 
(A7c) 
(A7d) 
(A7e) 
(A7f) 

The EJ21 scheme leads to a more general solution to the balance equations, and as 
in IW all the arbitrary parameters that appear in its solution can be associated with a 
residual choice of gauge. It has been explored in detail upto 4.5PN and the results are 
summarised below. We list the new general reactive solutions and the corresponding gauge 
transformations for the arbitrary parameters they contain. For brevity, the solutions are 
presented in the form : 'New solution' = 'Old solution' + 'Difference'. 

3. The 3.5PN and 4.5PN reactive solutions 

The reactive acceleration is assumed to have the following general form 



p' - 




AO, 


p' - 






P' - 

03 — 


0. 




Qoi = 






Q02 = 


0, 




Qo3 = 







1 111 

a = -^-[-(-42.5 + .A3.5 + A. 5 )rn + (£ 2 .5 + £3.5 + #2. 5 )v], (A8) 

V 



with ^4 2 .5 and £> 2 .5 given in Eqs. ( |A1| ) and ( |A4| ) and .A3. 5, B 3 _ 5 , A4.5 and £> 4 . 5 given by 
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With this form of the acceleration we have at 3.5PN 
dE* 8 mTm x ' * 15 



dt 

dJ* 
~dT 



rn 
15 V 



l2v 2 -llr 2 )+Y. n T ] yf ] 



i=i 



& ^ m 

5 r 2 



^(2^ 2 + 2^-3r 2 )+5: i S'r^., 



10 



,[3.5] ^[3] 



i=l 



(A9d) 

(AlOa) 
(AlOb) 



where y\ is given by Eqs. ( |2.15| ), 
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v , v — , v r , v — ,v — r , v r , — , — r, — r , r 



m 



.4-2 2 



/ m 



\ r 



,2 m -2 2-4 



/ m\ 3 /m 



V r 



V r 



„-,2 



(All) 



and TZ f , S'f'^ consist of corresponding linear combinations of the parameters involved. 
Repeating the procedure explained in the text, the 3.5PN reactive solution obtained is: 





= - 3t?)/% - 3^2 , 


(A12a) 




= h -\{2l + 39t7)/? 3 + 3p 2 , 


(A12b) 


fz 


= 2(1- 3 V )f3 3 + 4p 2 -5p 4 , 


(A12c) 


fi 


= fs + ^(56 + 15t7)/? 3 + 2p 2 + 5p 4 , 


(A12d) 


ti 


= 6p 4 , 


(A12e) 


fl 


= / 2 + (21 + 12^ 3 , 


(A12f) 


fl 


= / 5 -4r//3 3 , 


(A12g) 


fs 


= / 4 - 377/3 3 , 


(Al2h) 


ft 


= /6 ) 


(A12i) 


fio 


= o, 


(A12j) 


/ 

9i 


= 0, 


(A12k) 




= 9i , 


(A121) 


9', 




(A12m) 
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^ = 53-^(21 + 33^3 + 3^2, 
g' 5 = 2(l-3r ] )p 3 + Ap 2 -5p i , 

96 = 92, 

97 = 95 + ^(56 + r))[3 3 + 2p 2 + 5p 4 
g' 8 = g 4 + {21 + 9r})(3 3 , 

99 = 96, 

g'w = 6 P4 , 



(Al2n) 

(A12o) 
(A12p) 

(A12q) 

(A12r) 
(Al2s) 
(Al2t) 



where fi, gi are given by Eqs. (|2.9|) . The solution corresponding to Eqs. ( |2 . 1 0|) remains 
identical. 

Similarly at 4.5PN we have 



dE* 



dJ* 
~dT 



m 
15 ' r 



12v 2 — llf 2 I + 



28 
m 



(785 - 852//V 



+2(-1487 + 1392r])f V + 160(— 17 + t])—v 2 + 3(687 - 620r/)r 



r 



+8(367 - I577)— r 2 + 16(1 - Arj) — 



15 



i=i 



8 £ m 

5 r 2 



-.-,2 













.84 



(Al3a) 



2v + 2 r 

4(58 + 95t])—v 2 + 3(95 - 360r/)r 4 + 2(372 + 197?/)— r 2 



2-2 



(307 - 548r?K + 6 (-74 + 277r/) u 2 r 



+2(-745 + 277) - 



15 



+E#- 51 ^ [4] 



i=i 



(A13b) 



where 

,[5] 



3f J (i = 1...21) 



8 O-Z 4 (_) ^(_) ^ 



in 



.8-2 „,6 



/ m 



V r 



2 „.6-4 4 



/ m 



4^.4 4 .g ; 

f — r ,v r , v 



,m 



2j,8 



v — r , v r 
r 



m 



>v 

m 
r 



r J 



1 -2 2 

r ,v 



- r y 



r / 



V r 



8 ^10 



r , — r , r 
r 



(A14) 



and yf ] is given by Eq. Q2.15Q . 

Here ,iS'| 4 ' 5 ^ consist of linear combinations of the parameters involved. The 4.5PN 

reactive solution reads as: 



K = - ^(3 - 27r/ + 637/ 2 )/3 3 + ^(1 - 3r/)p 2 - 3 X 2 , 

o / 

h' 2 = ~(1 - 977 + 2l77 2 )/3 3 - 2(1 - 3 V )p 2 + 1(1- 3r / )p 4 + 4 X2 - 5 X4 
&s = hx + ~(3 - 207r/ + 7577 2 )/3 3 + ^(21 + 39r/)p 2 + 3 X2 , 



(Al5a) 
(Al5b) 
(A15c) 
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h' 4 = -3(1 - 3r])p 4 + 6x4 - 7 X 7 , (Al5d 

ti 5 = h 3 + (18 + 967] + 18r] 2 )(3 3 - (21 + 12ry)p 2 , (Al5e 

h' 6 = h 2 - 1(24 - 3977? + 9577 2 )/3 3 - 1(70 - llr l )p 2 + 1(35 + 657/)p 4 + 4 X2 + 5 X 4 , (A15f 

/4 = 8 X 7, (A15& 

K = h A + ~ (-353 + 195r7)77/3 3 + r/p 2 - J (84 + 25r/)p 4 + 2 X 4 + 7 X 7 , (Al5h 

o Z 

ti 9 = h 6 - 1(260 + 11977 + 30r/ 2 )/? 3 - (14 + 5r/)p 2 - (35 + 20r/)p 4 , (Al5i 

&io =h s - 1(306 + 489r/ + 48r/ 2 )/? 3 , (Al5j; 

h' n = h 10 , (A15k; 

/i' 12 = /i 9 - 1(12 + 87?7 - 24r7 2 )/? 3 , (A151 

ti l3 = h 8 - 1(8 + 49r/ + 34t7 2 )/3 3 + 2t/P2 + 5 V p 4 , (A15m 

/i' 14 = h 7 + 6t/[(1 - 3 V )(3 3 + p 4 ] , (Al5n 

h' lh = , (A15o 

k[ = , (A15p 

fc 2 = -1(3 - 277/ + 63t/ 2 )/? 3 + ^(1 - 37/)p 2 - 3 X 2 , (A15 q ; 

fea = ifei , (Al5r 

fc 4 = 1(1 - 97/ + 21r/ 2 )/3 3 - 2(1 - 3r/)p 2 + ^(1 - 377)p 4 + 4 X2 - 5 X 4 , (A15 S ; 

K = k 3 , (A15t 

k' 6 = k 2 + |(1 - 8I77 + 13t7 2 )/? 3 + 1(21 + 337/)p 2 + 3 X 2 , (Al5u 

o / 

k' 7 = -3(1 - 3r/)p 4 + 6x4 - 7X7 , (Al5v 
k' 8 = h - 1 (24 - 42177 - V 2 )fo ~ 1 (70 - 2hr 1 )p 2 + 1(35 + 55r;)p4 + 4 X2 + 5 X 4 , (A15 W ; 

K = k e + \( M + 525 ^ + 54t7 2 )/? 3 - (21 + 9?7)p 2 , (Al5x 

fc' 10 = k 5 , (Al5y; 

k'n = kw , (Al5z 

k' 12 = k 9 - - (159 + 28877 + 12r7 2 )/5 3 , (Al5aa 

k[ 3 = k 8 - 1(144 + 17977 + 40t7 2 )/3 3 - (14 + 67/)p 2 - (35 + 1577)p 4 , (A15bb; 

k' u = k 7 - 1(317 + 10577)77/33 - 1(84 + 377)p 4 - 3r/p 2 + 2 X4 + 7 X7 , (A15cc 

^5 = 8X7, (A15dd; 

where h iy ki are given by Eqs. (|2.16| ) of the text and Eqs. ( |2.17D remain the same. 
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4. The 3.5PN and the 4.5PN gauge arbitrariness 



Finally it can be shown that all the arbitrary parameters in the reactive solution may be 
absorbed in a choice of 'gauge' of the form 



m 



5x = -rj— (f' 5 + /a.5 + fi 5 )rx + (g' 2 _ 5 + g' 35 + si 5 )rv 



(A16) 



where f 25 and g 25 are given by Eqs. (|A6| 



/3.5 
93.5 
f±.5 



94.5 



P> 4 + P 9 > 2 - + Pf>,v 2 r 2 + Pt, 



yi „.2 
22 



m 



I ,.2^2 
23' 



') while fs.B> /4.5> #3.5 and #4.5 haVe the form 
2 



24" 



P' 



25 



+ ^ 
2 



Q> 4 + fe 2 - + Q' 23 ^ 2 + Q^-r 2 + q' 2 



25 



r 



r 



^4> 2 r 4 + Pi 7 -r 4 + PI 8 (J 



47 



18 



r 2 



2 



Q 41 ^ 6 + Q'v-v A + Q'^r 2 + Q' u v 2 (-) + Q' i5 
r \ r J 



-< 4 10' 



/ 2 ™ • 2 

v —r 



m 



+Q> r 4 + Q 47 -r 4 + Q' 4S 



r 2 + Q' 4 



19 



(A17) 



At 3.5PN we have 



P' 

r 2l 



P' 
r 22 



P' 

r 23 



P' 

r 2i 



P' 

r 25 

P' 
Q2I 
Q22 

Q23 

Q24 

Q25 

Q26 



P21 + ^(3 - 10t7)/3 3 + 7^(20p 2 + 17p 4 ) , 



-jP4, 

^23 + ^(5r//3 3 + P4), 

j P22 + ^(2 + 25r / )/5 3 -^(p 2 + P4) 

0, 
0. 

Q21 
1 
4 

Q23 



I(l_ 3r?)/g3 _ l(10 P2 + 7p 4 ), 

(1- 3rj)p 3 + i(2p 2 + p 4 ), 

i( 3 _ 877 ) /3 3_l(10p 2 + i3p 4 ) 
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Q22- 
1 

7^4- 



^(2 + 25r ? )/3 3 + i(p2 + p 4 ) 



(Al8a 
(A18b 
(A18c 
(A18d 

(Al8e 

(A18f 
(A18g 

(A18h 
(A18i 
(Al8j 
(A18k 
(A181 
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Similarly at 4.5PN we have 

PL = ~(1 - 977 + 21 V 2 )Ps + \{l ~ 3t7)P2 + ^(1 - 3r?)p 4 

-^(12x2 + 6x4 + 6x7), (A19a) 

PL = - ^( 1155 " 4817? ? + 36 V)& - 1(130 - 318 V )p 2 - 1(53 - 11777) P 4 

+ t^(105x2 + 84x4 + 68x7) , (A19b) 
105 

PL = \(1- 3^)P4 - Ife + 4 X7 ) , (A19c) 
PL = P42 + ^(420 - 1917^ + 967t7 2 )/3 3 + 1(55 - 22877)p 2 + 1^(220 - 834r7)p 4 

-^(15X2 + 14X4 + 13X7), (A19d) 
PL = P43 - ^(47 + 48^)^ - \vP2 ~ 1(8 - 17t ? ) P4 + yl(21X4 + 32 X7 ) , (Al9e) 

PL = -\xr, (A19f) 

27 11 
^7 = ^46 " gg^ 1 " 3 ^3 - - a VPa + -X7 , (A19g) 

P 48 = A 5 + -1(300 + 49357^ - 1360t7 2 )/5? 3 + ^(12^ - p 4 ) - 1( X4 + 2 X? ) , (Al9h) 
oOO zO 15 

PL = P44 - 1(92 + 12177 + 309t7 2 )/? 3 + 1(1 + 52?7)(p 2 + p 4 ) + ^( X2 + X4 + Xr) , (Al9i) 

^4io = 0, (A19j) 
Q 41 = , (A19k) 

Q42 = Q41 + 7^77(105 - 65977 - 347t7 2 )/? 3 + ^(1 - 377)p 2 + 1(1 - 377)p 4 
1 19 

-7373(10x2 + 7x4)-— X7, (A191) 
QL = ^(1 " 977 + 2l77 2 )/? 3 - ^(1 - 377)(2p 2 + p 4 ) + 1(12x2 + 6x4 + 4 X7 ) , (A19m) 
Ql 4 = Q42 ~ 7^(420 - I6O477 + 1434t7 2 )/3 3 - 1(80 - 30l77)p 2 - 1(40 - 13l77)p 4 

+ 1(10X2 + 9X4 + 8X7), (Al9n) 
QL = Q43 + 1^(175 - 63977 + 146t7 2 )/3 3 + 1(50 - 9977)p 2 + 1(94 - 18377)p 4 

-1(14x 2 + 14x4 + 12x7), (Al9o) 
QL = -^(1 " 3t7)P4 + 1(6X4 + 4 X7 ) , (A19p) 
QL = Q 46 + t!^7(121 - 363 V )[3 3 + ^ V p 2 + 1(5 - 877)p 4 - 1 %4 - , (A19q) 
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Q 48 = Q45 - 1(135 - 6477 - 11t7)/3 3 - 1(30 - 119r/)p 2 - 1(15 - 79r/)p 4 

+ 1(5X2 + 6X4 + 7 X7 ), (A19r) 

Q49 = Q44 + 1(92 + 12177 + 309r/ 2 )/5 3 - 1(1 + 52r/)(p 2 + p 4 ) - |( X2 + X4 + Xr) , (Al9s) 

Q' m = \x7- (Al9t) 



In the above, the P a b and Q a & are given by Eqs. ( |4.3| ) and ( |4.6| ) of the text. 

To conclude: the far-zone flux formulas and the balance equations by themselves do not 
constrain the reactive acceleration to be a power series in 777,1 and m 2 , or equivalently non- 
linear in the total mass m, as assumed in the paper, following IW. They are also consistent 
with the more general form of the reactive acceleration discussed in this appendix. 
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TABLES 



TABLE I. Comparison of four Alternative Schemes : IW21, IW22 (Minimal) IW11 and IWOO. 
N denotes the order of approximation, NV the number of variables, NC the number of constraints 
coming from balance equations, ND the number of degenerate equations, NI the number of inde- 
pendent equations and NA the number of arbitrary parameters determining the solution. In the 
NV column, a + b + c means a variables of reactive acceleration, b in energy ambiguity and c in 
angular momentum ambiguity. 



N 


NV 


NC 


ND 


NI 


NA 


IW21: IW Scheme 


2.5PN 


6+3+3 


12 


2 


10 


2 


3.5PN 


12+6+6 


20 


2 


18 


6 


4.5PN 


20+10+10 


30 


2 


28 


12 


IW22: Minimal Scheme 


2.5PN 


6+3 +1 


9 


1 


8 


2 


3.5PN 


12+6+3 


16 


1 


15 


6 


4.5PN 


20+10+6 


25 


1 


24 


12 


IW11 Scheme 


2.5PN 


6+6+3 


16 


3 


13 


2 


3.5PN 


12+10+6 


25 


3 


22 


6 


4.5PN 


20+15+10 


36 


3 


33 


12 


IWOO Scheme 


2.5PN 


6+10+6 


25 


5 


20 


2 


TABLE II. Comparison of the four Alternative Schemes: EJ21, EJ22, EJ11, and EJOO at 2.5PN 


level. The notation is as in Table [ 


|. In the NC column, a + b indicates that a constraints 


arise 


from enerj 


U balance and b from anj 


$ular momentum balance. 








Scheme 


NV 


NC 


ND 


NI 


NA 


EJ22 


12+3+1 


10+6 


2 


14 


2 


EJ21 


12+3 +3 


10+6 


1 


15 


3 


EJ11 


12+6+3 


10+6 


1 


15 


6 


EJOO 


12+10+6 


15+10 


3 


22 


6 
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